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ABSTRACT 

We study the location of massive disk galaxies on the Tully-Fisher relation. Using a 
combination of K-band photometry and high-quality rotation curves, we show that 
in traditional formulations of the TF relation (using the width of the global Hi pro- 
file or the maximum rotation velocity), galaxies with rotation velocities larger than 
200 km s -1 lie systematically to the right of the relation defined by less massive sys- 
tems, causing a characteristic 'kink' in the relations. Massive, early-type disk galaxies 
in particular have a large offset, up to 1.5 magnitudes, from the main relation defined 
by less massive and later-type spirals. 

The presence of a change in slope at the high-mass end of the Tully-Fisher relation has 
important consequences for the use of the Tully-Fisher relation as a tool for estimat- 
ing distances to galaxies or for probing galaxy evolution. In particular, the luminosity 
evolution of massive galaxies since z ~ 1 may have been significantly larger than es- 
timated in several recent studies. 

We also show that many of the galaxies with the largest offsets have declining rotation 
curves and that the change in slope largely disappears when we use the asymptotic 
rotation velocity as kinematic parameter. The remaining deviations from linearity can 
be removed when we simultaneously use the total baryonic mass (stars + gas) instead 
of the optical or near-infrared luminosity. Our results strengthen the view that the 
Tully-Fisher relation fundamentally links the mass of dark matter haloes with the 
total baryonic mass embedded in them. 

Key words: galaxies: spiral - galaxies: elliptical and lenticular, cD - galaxies: fun- 
damental parameters - galaxies: kinematics and dynamics - galaxies: statistics 



1 INTRODUCTION 



The notion of a tight correlation between absolute lumi- 
nosities of spiral galaxies and their rotation al velocities has 
been with us for about thirty years now l|Tullv fe Fisherl 
1 19771 ). This 'Tully-Fisher' (TF) relation has been con- 
firmed to hold over many decades in luminos ity (|Courteaul 



Il997l : iMcGaugh efail l2000l ; IVerheiienl l200ll) and in dif- 

at all ! 



and 



feren t galaxy environments (|Giovanelli et al.1 1 19971 : IWillickl 
1999). Barred an d unbarred galaxies s hare the same TF 
relation as well l|Courteau et alj |2003|). as do hig h 
low surface brigh tness galaxies ( Spravberrv et all 1 19951 ; 
IZwaan et al.l Il995l ) . Although spirals of diff erent morpho 
logical types do fol l ow different TF re l ations llRobertsI 197S 



2002 



Rubin et all 1 19851 : iHinz et al.l l200ll : iMathieu et all 
Russelll 12004 ). these offsets disappear almost entirely when 
using n ear-infrared photometry, rather than optical lum i- 
nosities (|Aaronson fc Mouldll 19851 : iPeletier fc Wirinerlll993h . 
indicating that most of the differences at bluer wavelengths 
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can be attributed to variations in star formation history 
along the Hubble sequence and the resulting differences in 
stellar populations. Similarly, when dynamical models are 
used to derive circular velocities for elliptical galaxies, it is 
found that the latter follow a similar TF relation to spirals, 
but offset to lower luminosities. This offset can again be ex- 
plained from differences in the stellar ma ss-to-light ratios 
jGerhard et al.ll200ll : Ide Riicke et~alll2007l ). 

The Tully-Fisher relation has become one of the most 
widely used relations in extragalactic astronomy. It has been 
commonly us ed as a powerful tool to estimate distances to 
galaxies fe.g. iTullv fc Pierce! 120001 and references therein). 
As a statistical correlation between fundamental properties 
of spiral galaxies, it has also been used to constrain nu- 
meric a l simulations of galax y form at ion ( Dalcanton et al.l 
19971: iNavarro fc Steinmet j l2000l: iBullock et al.l l200lh 
Pizagno et alj 120061 : iPortinari fc Sommer-Larsenl 120071 ). 
robe galaxy evolution on cosmological timescales 



Vogt et al l [1997: Ziegler et al. 2002; Milvang- Jense n et al 
20031 : [Bohm et alj|2004 iBamford et al.ll2006l : IWeiner et al. 
2006; Ichiu et alj|2007l ). or study the structure and stellar 
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populations of nearby galaxies dFranx fc, de Zeeuwl Il992l; 
ICourteau fe Rixll 19991 ; iBell fe de Jondl200ll ). 

Traditionally, the Tully-Fisher relation was studied us- 
ing the width W of a galaxy's global neutr al hydrogen line 
profile as a probe for its rotational velocity |Tullv fe Fisherl 
119771 ). With larger samples of spatially resolved kinematics 
(rotation curves) becoming available, the rotation velocities 
can be directly measured and recent studies tend to use the 
maximum rotation velocity (7 max ) instead. This subtle shift 
has the added advantage that the TF relation can be mea- 
sured at cosmological distances, using optical spectroscopic 
observations of high-redshift galaxies (see references above) . 

Most spiral galaxies have flat rotation curves and for 
these systems, both kinematic measurements (W and V ma x) 
generally yield nearly equivalent velocities. However, this 
is not the case for all galaxies. In galaxies whose rota- 
tion curves are rising or declining at large radii, it is not 
clear a priori which k inematic parameter yields the best re- 
sults. Vcrhciicn (2001, hereafter V01) presented an extensive 
study, based on a sample of 31 galaxies in the Ursa-Major 
cluster with well-defined rotation curves and K-band pho- 
tometry, of the influence of the shape of a rotation curve on 
the location of a galaxy on the Tully-Fisher relation. Galax- 
ies with declining rotation curves were found to lie system- 
atically on the high-velocity side, when using the maximum 
rotation velocity Vm ax or the width of the Hi profile. Since 
these are typically bright systems with high rotation veloci- 
ties, this result suggests that there may be a change in slope 
in the TF relatio n at the high-luminosity end, a possibility 
also hinted at bv lPeletier fc Willnerl (|l993l ). Clearly, the ex- 
istence of a 'kink' in the Tully-Fisher relation has important 
consequences and, if not corrected for, will lead to system- 
atic biases when deriving cosmological distances or probing 
galaxy evolut ion. 

However, IVOll also showed that the systematic offset at 
the high mass end disappeared when, instead of Vmax, the 
asymptotic rotation velocity Vksymp in the outer, flat parts 
of the rotation curve was used. Bearing in mind that the 
asymptotic rot ation velocities are determined by the dark 
matter haloes, IVOll interpreted this result as a strong indi- 
cation that the TF relation fundamentally links the total 
baryonic content of (disk) galaxies with their dark matter 
haloes, and that as far as this relation is concerned, there is 
no fundamental difference between low and high mass galax- 
ies. 

Unfortunately, the sample of lVOll only contained a small 
number of galaxies with declining rotation curves and the 
declines were modest. Thus, the change in slope of the high- 
luminosity end of the TF relation when going from V ma x 
to \4symp was small and his results only gave a first indi- 
cation for how to interpret the curvature at the high-mass 
end of the TF relation in terms of the kinematics of the 
galaxies. Recently, however, a sample of high-quality rota- 
tion curves for high-luminosi ty early-type disk galaxies (SO - 
Sab) has become available l|Noordermeer et al.ll2007l , here- 
after N07). This sample contains many strongly declining 
rotation curves where the difference between the maximum 
and asymptotic rotation velocities is much larger than in 
Verheijen's galaxies. In this paper therefore, we combine 
Noordermeer's sample with the data for th e galaxies with 
flat and declining rotation curves from I VP ll . with the aim of 
further investigating the change in slope at the bright end 



of the TF relation and its relation with the shapes of rota- 
tion curves. We also include the recent s ample of rotation 
curves from ISpekkens fc Giovanellil (|2006l . hereafter S06) in 
our analysis; these are all massive late-type spiral galaxies 
with flat rotation curves and provide a valuable benchmark 
at the high-luminosity end of the TF relation. 

The outline of the remainder of this paper is as fol- 
lows. In section [2] we describe the selection of galaxies from 
the three different samples for the present study, together 
with the photometric and kinematic ingredients for the TF 
analysis. In section [3j we present the Tully-Fisher relations 
and show that there is indeed a change in slope around 
200 km s -1 in the relations with W and Knax- We also show 
that this 'kink' is largely, but not entirely, removed when we 
use the asymptotic rotation velocity Vasymp- In section [3] we 
search for possible explanations for the observed change in 
slope and show that, when we consider the total baryonic 
mass (stars + gas) in our galaxies, an entirely linear rela- 
tion is recovered. Finally, we discuss some implications of 
our results in section [5] 



2 SAMPLE SELECTION AND 
OBSERVATIONAL DATA 

2.1 Sample 

We comb i ne d ata fro m the three separate studi es by 
IVerheiienl lioOll. V011. [Spekkens fc Giovanellil l|2006l , S06) 
and iNoordermeer et all (|2007l . N07). The main criteria for 
inclusion of galaxies from these sub-samples in the present 
study is that their rotation curves are of good quality and 
that the maximum and asymptotic rotation vel ocitie s are 
well denned. We thus select the 22 galaxies from IVOll with 
flat or declining rotation curves. For galaxies whose rotation 
curves rise until the last measured point, the maximum and 
asymptotic r otatio n velocities are ill-defined, so we excluded 
galaxies from I VP ll with such rotation curves. We include all 
galaxies from the sample of massive, late- type spiral galaxies 
from ISPrj . except UGC 2849, which appears to be interact- 
ing with a companion and whose kinematics are lopsided. 
The remaining 7 galaxies have well-defined rotation curves, 
some of which appear to be weakly declining in the outer 
regions. However, given the uncertainties in the inclination 
of the gas layer at large radii, they are also consistent with 
bein g entirely flat. Finally, we include all 19 galaxies from 
IN07I . These are all early-type disk galaxies (SO - Sab), and 
mostly massive (Vmax <^ 200 km s -1 ). Many of these sys- 
tems have declining rotation curves, with large differences 



Our final sample consists of 48 galaxies. They are listed, 
together with the photometric and kinematic data, in ta- 
bled] 



2.2 Photometric data 

We have chosen to study the Tully-Fisher relation using K- 
band photometry. Near-infrared photometry has a number 
of clear advantages over conventional optical data. Foremost, 
the effects of dust extinction is strongly reduced. Several 
galaxies in our sample are highly inclined with respect to the 
line of sight, such that a large fraction of the light is absorbed 



Massive galaxies on the Tully-Fisher re 



Table 1. Photometric and kinematic properties use d for the Tully-Fisher analysis: (1) galaxy name; (2) total 
apparent K-band magnitude and error, taken from T96 (Verheijen's galaxies) and the 2MASS galaxy cata- 
logue (Spekkens' an d Noordermeer's galaxies); (3) correction for galactic foreground extinction, taken from 
ISchlegel et al,l (l998j); (4) correction for internal extinction, calculated using equation [T] (5) assumed distance; 
(6) magnitude error due to distance uncertainty; (7) resulting absolute magnitude; (8) maximum and (9) asymp- 
totic rotation velocities from the rotation curve and (10) width of global Hi profile, corrected for instrumental 
broadening, random gas motions and inclination. 

galaxy 7^ If A*~ D <5A/ dist m| V max V asymp W^ R 
mag Mpc mag kms -1 



(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) 
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1~The rotation curve of UGC 11914 only extends out to about 3.3 R-band disk scale lengths and does not 
show the characteristic decline that is seen in other galaxies of similar type and luminosity (see N07) . Since 
in many comparable cases, the decline in the rotation velocities sets in outside the optical disk only, it is well 
possible that the rotation curve in UGC 11914 would also decline at larger radii if we were able to measure 
it. The asymptotic rotation velocity Vksymp for UGC 11914 seems therefore ill-defined and we excluded this 
galaxy from the sample for the asymptotic velocity relations. 
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Figure 1 . Com parison of 2MASS K-band magnitudes with the 
data from lT96l . Data points show the difference betw een th e total 
apparent magnitudes from 2MASS and those from T96; error- 
bars show the combined errors. Open symbols indicate low sur- 
face brightness galaxies (/irj,A' > 18 mag arcsec -2 ), filled symbols 
indicate HSB systems. Black and red vertical arrow s ind icate the 
2MASS apparent magnitudes for the galaxies from |N07| and ISOtj 
respectively. 



by internal dust. In the optical bands, this effect can be 
very significant (<; lmag), leading to large uncertainties in 
the derived absolute luminosities. In the near infrared, these 
uncertainties are less severe. In addition, recent bursts of star 
formation can have a large effect on the optical luminosity 
of a galaxy, whereas the K-band is much less sensitive to 
such events. 

Deep K- band photometry f or the galaxies from IVOl] is 
available from lTullv etafl (|l996l . hereafter T96). NeitherlSOd 
nor lN07l had private near-infrared photometry at their dis- 
posal, but all galaxies from these subsamples have been ob- 
served in K in the framework of the Two Micron All Sky 
SurvejE Although 2MASS is a very homogeneous survey, it 
is also less deep than the data from |T96| . To judge the reli- 
ability of the 2MASS photometric data, and to ensure that 
we do not introduce systematic differences between our var- 
ious subsamples, we compared the K-band magnitudes from 
T96 with the total extrapolated magnitudes from 2MASS, 
for all systems which were detected in both surveys. 

Figure [T] shows that there is good agreement between 
the 2MASS photometry and Tully's data for high surface 
brightness (HSB) galaxies. For the galaxies with measured 
central surface brightness fJ,o,K < 18.0 mag arcsec -2 , the av- 
erage difference between 2MASS' and T96's magn itude s is 
0.01 mag, with a standard deviation of 0.14 mag. Ivoil re- 
ported, based on an internal comparison of their observa- 
tions from different nights, an average photometric uncer- 
tainty of 0.08 mag in the K-band magnitudes of |T96| . This 
implies that the uncertainties in the 2MASS magnitudes 
must be approximately 0.11 mag, significantly larger than 



1 2MASS (http://www.ipac.caltech.edu/2mass/) is a joint 
project of the University of Massachusetts and the Infrared Pro- 
cessing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Administration 
and the National Science Foundation. 



the average errors given by the 2MASS database (typically 
0.03 mag for the HSB galaxies in figure [T}. Photometric 
errors of 0.03 mag seem somewhat optimistic, in particu- 
lar since the 2MASS images are significantly less deep than 
those of lT96l . 

Figure Q] also shows that for low surface brightness 
galaxies, the 2MASS magnitudes are not re liable and some- 
times deviate strongly from the values from lT96l . Inspection 
of the 2MASS images for these galaxies shows that they are 
of such low surface brightness that they are often barely de- 
tected (in fact, most of the LSB galaxies from Tully's sample 
are not detected at all by 2MASS). It appears that 2MASS 
has missed a large fraction of the fl ux of many LSB galaxies, 
compared to the deepe r dat a fro m |T9q . 

The galaxies from lN07l and ISOq are all luminous, high 
surface brightness galaxies. For these objects, it seems safe 
to adopt the 2MASS magnitudes, but we assume an average 
photometric uncertainty of 0.11 mag instead of the smaller 
errors given by 2MASS. 

To convert the apparent magnitudes to absolute lu- 
minosities, we first correct for extinction effects. Correc- 
tions for Galactic fo reground extinction were taken from 
ISchlegel et alj (|l998l ); they are generally small (< 0.1 mag 
for 85% of the galaxies). Several correction schemes exist to 
determine the amount of internal extinctio n; the most com- 
monly use d methods are ba s ed on those of iTullv fc Fouauel 
1 19851 ) or iGiovanelli et~ai1 l| 19941 ). Here, we follow IVOU 
who employed th e following relation, originally derived by 
ITullv et al. (|l998l ). for the internal extinction parameter: 



-7Alog(£), 



(1) 



with b/a the observed minor-to-major axis ratio of the op- 
tical image . 7a is wavelength dependent and was found by 
ITullv et alj (| 199ST ) to depend also on the absolute luminos- 
ity of the galaxy: brighter galaxies contain on average more 
dust than fainter ones. These authors used the Tully-Fisher 
relation itself, in an iterative way, to express the absolute lu- 
minosity in terms of the Hi line width, and give the following 
description for 7a: 



7a = ax + /^(logWa 



2.5), 



(2) 



where VKjo'h is the Hi line width, corrected for inclination 
and broadening due to instrumental effects and random gas 
motions (see below). For the wavelength depen dent param- 
eters cea and /3a, we use the values given by ITullv et al.l 
(1998). Note that even in the K-band used here, we find sig- 
nificant internal extinction corrections, although they do not 
exceed 0.1 magnitude except in the most inclined systems. 

Distanc es fo r our galaxies were derived as follows. All 
galaxies in IVOll belong to the Ursa Major cluster and 
have therefore roughly equal distances; throughout this 
study, we adopt Ver heijen's value of 18.6 M pc (ta ken from 
ITullv fc Pierce] I2000I ). For the galaxies from l N07l. w e used 
the distances as derived bv lNoordermeer et al. ( 20051 ) . which 
were based on a simple Hubble-flow model, corrected for 
Virgo-centric infall and assuming a H ubble constant of 
75kms -1 Mpc -1 . For the galaxies from Isod we used the 
same methocjf]. 



Note that our distance estimates for the galaxies from ISOa are 
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In table [T] we give for each galaxy the raw apparent 
magnitudes and photometric errors, the galactic foreground 
and internal extinction parameters, A J and A l K , the as- 
sumed distances and the resulting corrected absolute mag- 
nitudes used in the remainder of this paper. 

Finally, in addition to the photometric errors, we also 
account for the uncertainties in the absolute magnitudes 
which arise fr om e rrors in our distance estimates. For the 
galaxies from IVOll these are defined by the depth of the 
Ursa Major cluster; here, we adopt Verheijen's estimate 
of 0.17 mag. For the other galaxies, we estimate the dis- 
tance uncertainties by assuming a typical peculiar velocity 
of 200 km s" 1 . The resulting errors on the absolute magni- 
tude <5Mdist are given in column (7) of table [T] and were 
added quadratically to the photometric errors to obtain the 
total uncertainties used for the subsequent analysis. 

Note that we have not included the uncertainties in the 
internal extinction corrections in our error budget here. Es- 
pecially for the highly- inclined systems in our samples, the 
corrections A\ are significant and galaxy-to-galaxy varia- 
tions of 7a can cause substantial deviations from the cor- 
rections used here. However, without additional information 
about the dust content of our galaxies, the corrections listed 
in table [1] are the best possible estimate and the uncertain- 
ties are very difficult to quantify. We will ignore these un- 
certainties for the subsequent analysis here, but note that 
some of the scatter in our TF relations might be explained 
by the uncertainties in the internal extinction. 



random gas motions. IVerheiien fc Sancisil i|200ll ) present an 
extensive discussion about the suitable choice for Wt\ here 
we simply copy their preferred value of Wt = 22 km s _1 . The 
derived profile widths W^o r were corrected for inclination 
using the values derived in the original papers; for warped 
galaxies we use the inclination in the inner regions where 
most of the gas is concentrated. 

The two other kinematic parameters are the maximum 
rotation velocit y Kn ax and the asymptotic rotation velocity 
V^symp (Vfiat in VOID, which are both derived directly from 
the rotation curves. 

For the galaxies from I VOll , we copied the errors on the 
kinematic par amet ers without further modification. For the 
galaxies from IN07I . we used the inclination uncertainty de- 
rived in the original paper to estimate the errors on the 
profile widths. The errors on the two other kinematic pa- 
rameters were estimated by eye, based on the errors in the 
rotation curves. The latter include contributions from mea- 
surement errors, kinematic asymmetries and inclination un- 
certainties; for the more face-on galaxies, the l atter are usu- 
ally dominant. Finally, for the galaxies from ISOrl we used 
the uncertainty in the ellipticity of the optical images, given 
in the original paper, to estimate the inclination uncertainty, 
but adopted a minimum of 5° to account for the possibil- 
ity of undetected warps in the outer parts of the gas layer. 
The errors on Vmax and Vasymp were again estimated by eye, 
based on the errors in the rotation curves. 

All three kinematic parameters and corresponding er- 
rors for each galaxy are given in table [1] 



2.3 Kinematic data 



We use three kinematic parameters for our Tully-Fisher 
analysis. The first is the inclination corrected width of the 
Hi line profile. For this, we used the instrumental bro a denin g 
corr ected measurements from Verheiien fc Sancis 
ISOd and iNoordermeer et all l)2005l ) respectively. We ap- 
ply an additional correction for broadening of the profiles 
caused by random gas motions in the gas disks, using the 
prescriptions given by IVerheiien fc Sancisil <(200ll . based on 
iTullv fc Fououil 19851 ): 



(^2o,h) 2 = (VK 2 c ) 2 + W t 2 



2c 



■(*)' 



(3) 



with W20 the profile width corrected for instrumental broad- 
ening. Wc indicates the profile width where the transition 
from a Gaussian t o a boxy shape occurs; we assume here 
W c = 120 kms -1 |de Vaucouleurs et al.lll983h . In practice, 
most of our galaxies have profile widths larger than W c , 
such that equation[3]yields a linear subtraction of Wt, which 
gives the amount by which a profile is broadened due to the 



slightly different from those given in the original paper. Part of 
the discrepancy comes from the fact that we assume h = 0.75 
instead of 0.70. Further differences arise from the treatment of 
Virgo-centric streaming motions, but are small (< 0.1 mag in all 
cases). 



3 TULLY-FISHER RELATIONS 

In figure [5] we show the Tully-Fisher relations for the c om- 
bined sample of galaxies. Compared to the study in lVOli our 
combined sample greatly increases the number of luminous 
galaxies, with our most luminous system, UGC 2487, being 
1.2 magnitudes brighter than the most luminous object in 
Verheijen's sample. 

We fitted direct linear relations of the following form to 
the combined data sets: 



M C K = M 2 .6 + S (log W - 2.6) . 



(4) 



where is the extinction corrected absolute magnitude 



and W is either W£ R , 2V m 



or 2T4sym P . We calculated the 



zeropoint M2.6 at log W = 2.6, rather than 0, to avoid strong 
covariances between the two free parameters. The fits were 
done using a simple x 2 - m i mm i sa tion procedure, taking the 
errors in both directions into account. They are indicated 
with the magenta lines in figure [2] for compar ison, we also 
show (in black) the fits to the 'F'-sample from I VOll . i.e. the 
fits made to galaxies with flat rotation curves only. 

The results of the fits to the combined sample are sum- 
marised in tabled The scatter around each of the fits given 
in this table is a weighted rms scatter, with the weight for 
each data point given by Wi = (p 2 M A + S 2 of og w^) -1 - Here, 
(JM,i are the magnitude errors, a\ gw,i are the errors in the 
parameters on the z-axes and S are the fitted slopes of the 
relations. 

A number of interesting results can be recognised from 
figure [2] and table [2] First of all, our data strongly sug- 
gest a 'kink' in the Tully-Fisher relation: the TF relation 
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Figure 2. Tully-Fisher relations in the K-band, using the corrected widths of the global Hi pr ofiles (left) and the maximum (middle) 
and asymptotic rotation velocities (right). Black and green data points show galaxies fr om | V0ll with flat and declining rotation curv es 
respectively; blue data points show galaxies from |S06| and red points show galaxies from lN07i The black lines show the fits from IVOll to 
galaxies with flat rotation curves, while the magenta lines show the fits to the combined sample (see tabled- 



Table 2. Results from the least-x 2 fits to the K-band Tully-Fisher relations shown in figure [2] 



kinematic 


# of 


M 2 .6 


S 


scatter 


^rcd 


Q 


parameter 


points 














mag 




mag 






vy 20,R 


48 


-23.69 ± 0.03 


-8.65 ± 0.19 


0.38 


2.89 


3.3 ■ 10- 11 


Vmax 


48 


-23.64 ± 0.04 


-8.77 ± 0.22 


0.44 


2.98 


7.2 ■ 10~ 12 


^asymp 


47 


-23.95 ± 0.04 


-9.64 ± 0.25 


0.36 


1.76 


5.4 ■ 10~ 4 


Vasymp (UGC 3993 and 6787 excluded) 


45 


-23.97 ± 0.04 


-9.71 ± 0.26 


0.34 


1.54 


6.8 ■ 10- 3 



seems to become shallower above a rotation velocity of about 
200kms _1 (equivalent to M % » -23.75). This i s cons istent 
with the claim first made bv lPeletier fc Willnerl l|l993h . but 
the effect is much more clearly visible here than in their 
data. The kink is most apparent when using the width of 
the global profile (W 7 ^ 1 ^) or the maximum rotation veloc- 
ity V ma x as kinematic parameter. In the left an d mi ddle 
pane ls in the figure, almost all galaxies from the |S06| and 
|N07| samples lie to the right of the relation defined by the 
inter mediate-mass galaxies with flat rotation curves from 
IVOll . As a result, the fits to the combined samples have a 
shallower slope than those from IVerheiier] to the galaxies 
with flat rotation curves. In particular, many of the mas- 
sive, early- type disk galaxies from lN07l lie far away from the 
rela tion defined by the later- type spirals from IVerheiienl and 
|S06| ; when expressed in terms of a luminosity offset, their 
distance from the main relation can be as large as 2 magni- 
tudes. 

However, the turnover is greatly reduced when using the 
asymptotic rotation velocity Vasymp as kinematic parameter. 
Many of the galaxies in our sample with V m ax > 200 km s^ 1 , 
in particular the early-type disks from lN07l have strongly 
declining rotation curves. Galaxies such as UGC 4458, 4605, 
3546, which lie far to the right of the main relation in the 
left and middle panels, shift to the left when the (lower) 
asymptotic velocity is used, thereby straightening the TF 
relation. This effect is also reflected in the values for the rms 
scatter and \ 2 of the data around the fits and the 'goodness- 
of-fit' parameters Q. Table [2] shows that the TF relations 



using T4symp are much better represented by straight lines 
than the relations with the other two parameters. 

It is important, however, to note that the kink does 
not disappear completely when using Vasymp as kinematic 
parameter. Even in the right hand panel of figure [5] most 
galaxies at the bright end lie to the right of the relation 
defined by the l ess luminous galaxies with flat rotation 
curves from IVOll . indicating that the slope of the TF re- 
lation also changes when the asymptotic rotation velocities 
Vasymp are used. This is also confirmed when we fit linear 
relations to the subsamples of galaxies fainter and brighter 
than MJf ~ —24; in this case, we find that the slope changes 
significantly, from —10.5 ± 0.5 to —7.7 ± 0.9 respectively. 

The scatter in our relation using the asymptotic rota- 
tion velocities is sys tema tically larger than the correspond- 
ing value found by IVOll . This difference can partly be at- 
tributed to the fact that Verheijen's galaxies all lie in the 
Ursa-Major cluster, such that the distance uncertaintie s are 
sma ll (at least in a relative sense) . The galaxies from IS06l 
and Inot lie predominantly in the field, where peculiar mo- 
tions with respect to the Hubble flow lead to errors in the 
derived distances and absolute luminosities (see table [TJ, 
and thus to additional scatter in the TF relations. More im- 
portantly however, the scatter in our relations is artificially 
increased by the deviations from a straight line. Even in the 
relations using the asymptotic rotation velocity Vasymp, the 
kink around 200kms~ 1 cause s sys tematic deviations from 
the fitted straight lines. Since IvOl] had only 2 galaxies with 
Vksymp > 200 km s , the change in slope at the high mass 
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Figure 3. Deviations of our objects with respect to the K-band 
vs. Vksymp TF relation from table[2] Symbol colours are the same 
as in figure [2] Errorbars are the effective errors, calculated by 
combining the magnitude and velocity uncertainties (see text). 
Dotted lines indicate the formal 3cr limits. 
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Figure 4. Deviations of our objects with respect to Verheijen's ri- 
band vs. Vksymp TF relation for galaxies with flat rotation curves, 
as a function of inclination. Symbol colours are the same as in 
figure [2] The dotted lines show the adopted internal extinction 
corrrections for fiducial galaxies with Hi profile widths of 200, 400 
and 600 km s -1 (top to bottom respectively). 



end had virtually no influence on the scatter around the fits 
for his data. 

Finally, the scatter in the relations with Vasymp is heav- 
ily influenced by two galaxies with unusually large devia- 
tions from the main relation: UGC 3993 and UGC 6787 
(3.4 and 3. Oct outliers respectively). In figure [3] we show 
the deviations of the galaxies from our sample with re- 



spect to the Mjf vs. Vis 



fit described in table [2] The 



errorbars in the figure take the uncertainties of the points 
in both directions into account, and were calculated as 
(A(M-M fl t)) 2 = (AM) 2 + (Alog(2V)-S') 2 ,withS= -9.64 
the slope of the fit. 

UGC 3993 is a nearly face-on galaxy (i ps 20°, see lN07T ). 
with correspondingly large uncertainties in the derived ro- 
tational velocities (also reflected in the large errorbars in 
figures [2] and [3]) . The deviation of this galaxy may be due 
partly to a slight under-estimation of the true inclination 
angle. It seems unlikely, however, that the inclination un- 
certainty is responsible for the full deviation of UGC 3993; 
to bring this galaxy to the middle of the relation would re- 
quire an inclination angle of 30° , a value which seems to fall 
outs ide the range supported by the observational data (see 

103). 

The most likely explanation for the offset of UGC 6787 
from the main relation lies in the distance determination. 
UGC 6787 lies close to the center of the Ursa Major cluster, 
but has a redshift that is about 200 k m s" 1 higher than the 
high velocity envelope of the cluster (|Tullv et al.lll996T ). It 
may thus lie behind the main cluster, which might in turn 
imply that it is being drawn into the cluster and that the 
recession velocity is lower than expected in the case of pure 
Hubble flow. In that case, our adopted distance of 18.9 Mpc 
and the derived luminosities are too small, explaining the 
offset in the TF relations. 

For completeness, we also list in table[5]the results of fits 
to the TF relations with the asymptotic rotation velocities 
when UGC 3993 and 6787 are excluded from the sample. 



It is clear that the exclusion of the two points from the fits 
does not lead to significantly different slopes or zeropoints. 
As expected, the scatter is reduced, but since the kink in 
the relations has not been removed by the exclusion of the 
two discrepant points, t he va lues are still larger than the 
corresponding ones from IVOll and the \ 2 an< i Q-parameter 
indicate that the deviations from a straight relation are still 
real. 



4 DISCUSSION 

4.1 Systematic variations along the TF relation 

What could be the origin of the change of slope in the Tully- 
Fisher relation? Does it truly indicate a break in the relation 
between baryons and dark matter, or could it be explained 
by other effects? An obvious possibility is that some of the 
corrections from sections 12.21 and 12.31 are inaccurate for the 
brightest galaxies. Such an effect may lead to a systematic 
deviation of these galaxies from the TF relation for lower- 
mass galaxies, precisely as observed. 



4-1.1 internal extinction 

The largest uncertainty in the corrections which we applied 
to our raw data lies, most probably, in the treatment of 
the inter nal extinction. We corrected for internal extinction 
following iTullv et al.l ljl998h . They derived their corrections 
from a sample of galaxies which spans a roughly similar 
luminosity range as ours and, importantly, observed no sys- 
temat ic deviations of t he br ightest galaxies. On the other 
hand, iGiovanelli et al.l fl995) did find a hint for excess ex- 
tinction in the brightest spirals. Given the uncertainties, it 
is useful to explicitly verify whether the deviations of our 
most massive galaxies with respect to the TF relation for 



8 E. Noordermeer & M. A. W. Verheijen 



1.5 



0.5 



- 


1 

H 

T 


"-ft 


I _ 
- 



















-20 



-32 



-24 



-26 



Figure 5. Colour-magnitude relation for all galaxies in our sam- 
ple for which B-R colours are available. Colour information for 
galaxies from | V0ll were taken from the same paper, whereas 
iNoordermeer &; van der Hu lst (2007) give colours for a subset of 
the galaxies from lN07l . The black line shows a least-squares linear 
fit to the data points. 



low- and intermediate-mass galaxies can be explained as be- 
ing due to an underestimation of the internal extinction in 
these systems. This is done in figure [4] where we plot the 
deviations of all galaxies with respect to IVerheiie'rH 's TF re- 
lation for spirals with flat rotation curves, as a function of 
inclination. 

Figure U shows a tentative trend for the galaxies from 
ISOrj (blue points), with the edge-on systems deviating more 
than the galaxies with intermediate inclinations. To explain 
this trend, about a factor three more extinction is required 
in th ese g alaxies than assumed in section 12.21 The galaxies 
from IN07I (red points), however, do not follow an obvious 
trend and strongly deviating cases are found among both 
near-face-on and highly-inclined galaxies. 

Based on our data, we cannot strictly rule out the possi- 
bility that the internal extinction in the most massive spiral 
galaxies is higher than we assumed. However, figure [4] shows 
that if this is indeed the case, it still cannot explain the 
deviations of the low-inclination galaxies. It therefore seems 
unlikely that the change of slope at the high-mass end of the 
Tully-Fisher relation can be explained purely as an effect of 
excess extinction in the most massive disk galaxies. 



4-1.2 stellar population variations 

Another possibility which could explain a change of slope 
in the Tully-Fisher relation is that the most massive galax- 
ies have stellar populations with higher mass-to-light ratios 
than the less luminous systems, such that they are underlu- 
minous for a given rotation velocity. To first order, such an 
effect is likely to be present in disk galaxies, given the well- 
known fact that bright galax ies tend to be redd e r than their 
less-l uminous counterparts l|Visvanathanlll98ll ; iTullv et al.l 
1982). This colour-magnitude relation is also present in our 
own sample: in figure [5] we plot the B-R colours as a func- 
tion of absolute luminosity, for the subset of our galaxies for 
which the former are available. 

The colour-magnitude relation itself cannot, however, 
be responsible for the kink in the TF relation. Figure[5]shows 



that the colour changes linearly with absolute magnitude. 
Since mass-to-light ra tio changes exponentially with colour 
(jBell fe de Jondl200ll ). this gradient will translate in a lin- 
ear effect in the Tully-Fisher relation. In other words, a lin- 
ear colour- magnitude relation will introduce a change in the 
global slope of the Tully-Fisher relation. Specifically, using 
the fitted slope of the colour-magnitude relation in figure [5] 
in conjunction with the slope of the colour vs. mass-to-light 
ratio relation from iBell fc de Jond . it can be derived that 
the stellar mass Tully-Fisher relation will be 10% steeper 
than the relations shown in figure [2] and table [2] Impor- 
tantly, however, a linear colour-magnitude relation cannot 
introduce a turnover in the TF relation. Only a break in 
the colour-magnitude relation at the same luminosity as the 
kink in the TF relation can explain the latter, but figure [S] 
shows that such a feature is not present in our sample. 



4.2 The Baryonic Tully-Fisher relation 

Having established that the observed change of slope in 
our K-band vs. V aS ym P TF relation is not due to varia- 
tions in internal extinction or the properties of the stel- 
lar populations, we now consider whether trends in rela- 
tive gas content play a role. In this context, it is inter- 
esting to consider the so- called 'Baryonic Tully-F isher re- 
lation', first dis cussed by iMcGaugh et al.l l|200u1 . see also 
iMcGaughl [20051 ). They showed that there exists another 
break in the Tully-Fisher relation at the low luminosity end 
(around Kot = 90kms -1 ), below which galaxies are also 
under-luminous. They were, however, able to restore a lin- 
ear TF relation when, instead of using the stellar luminosity, 
they adopted the total observed baryonic mass (stars and 
gas). Since dwarf galaxies contain on average more gas than 
high er-luminosity spirals ( relati v e to the optical luminosity , 
e.g. lHavnes fc Gioyaneilil Il984l; Verhciicn fc Sancisil 2001 
ISwaters et al.l |2002| ; iGeha et al l 120061 ). the former shifted 
upwards more than the latter, and the kink in the Tully- 
Fisher relation disappeared. Luminous galaxies, on the other 
hand, generally con t ain re latively little gas. In particular, in 
INoordermeer et al.l (|2005T ). i t was shown that some of the 
early- type disk galaxies from lN07l are very gas-poor. In such 
galaxies, the baryonic budget is dominated by the stars and 
adding the gas contribution will not lead to a significant 
increase in brightness. 

We have investigated whether the observed break at the 
high mass end in our Tully-Fisher relations can be explained 
as a result of the relatively low gas-content of the most lu- 
minous galaxies in our sample. To this purpose, we have 
converted the total K-band luminosities to stellar masses, 
assuming an average mass-to-light ratio of the stellar popu- 
lati ons of 0.8. The choice o f M*/ Lk = 0.8 was made follow- 
ing IMcGaugh et al.l l|2000h , and is consistent wi th the av- 
erage valu es found from maximum-di sk fits from IVerheiienl 
(|l997f ) and lPalunas fc Williams] (|2000T ). In reality, the values 
of M*/Lk are expected to vary from galaxy to galaxy, but 
the variations will be modest in the K-band and we ignore 
this variation for the moment. The possibility that the as- 
sumption of constant M/L gives rise to additional scatter in 
the baryonic TF relation is investigated below. 

For the gas masses, we used the to tal Hi ma sses as 
given by I VP J INoordermeer et"aL (2005) and ISOfj respec- 
tively, multiplied by 1.43 to account for the presence of he- 
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Figure 6. Baryonic Tully-Fisher relation. The total K-band luminosity was converted to stellar mass, assuming a mass-to-light ratio of 
the stars of 0.8 (see text), and added to the neutral gas mass. Colours of symbols and lines are the same as in figure[2] 



Table 3. Results from the least-x fits to the Baryonic Tully-Fisher relations shown in figure [< 



kinematic 


(logAf) a .6 


S 


scatter 




Q 


parameter 






mag™ 






yy 20,R 


10.78 ± 0.01 


3.04 ± 0.08 


0.32 


2.17 


2.2 • TO -6 


Vmax 


10.77 ± 0.01 


3.05 ± 0.09 


0.40 


2.57 


2.4 ■ 10" 9 


Vasymp 


10.88 ± 0.02 


3.36 ± 0.10 


0.36 


1.92 


7.7- 10~ 5 


V asymp (UGC 3993 and 6787 excluded) 


10.89 ± 0.02 


3.38 ± 0.10 


0.34 


1.72 


1.0 • 10" 3 


a For consistency with table |2| we express 


the scatter in 


magnitudes. 


The scatter 


in log(M bar ), as 



shown in figure [6] is a factor 2.5 lower. 



Hum. The stellar and gas masses were then added to obtain 
an estimate for the total baryonic mass content in our galax- 
ies. The resulting baryonic TF relations are shown in figure[6] 
for the same three kinematic parameters as used before; the 
resulting least-x 2 fits are summarized in table [3] 

As expected from the general trends of gas content with 
total luminosity mentioned above, the baryonic Tully-Fisher 
relations have a shallower slope than the standard stellar 
TF relations. But more importantly, the 'kink' in the TF 
relations from figure [2] is reduced. The relations with W^q r 
and Vmax are better represented by a straight line than the 
original relations in figure [2] the scatter and x 2 are reduced 
with respect to the original values and the Q-parameters are 
higher. Thus, our data show that the concept of the Bary- 
onic Tully-Fisher relation is not only useful to increase the 
linearity in the TF relation at the low luminosity end, but 
that in addition, it also reduces the kink around 200 kms -1 . 
However, this effect is smaller than the one discussed in the 
previous section, and the Baryonic Tully Fisher relations us- 
ing W£n R or Vmax are still worse than our original, stellar 
luminosity vs. Vasymp relation shown in figure[2]and table[2] 

The inclusion of the gas does, however, remove the small 
kink that was still present in the latter relation as well, and 
the baryonic mass vs. Vasymp TF relation appears to be con- 
sistent with a linear relation over the full extent of our data. 
However, at the same time, the scatter around the mean 
relation does not appear to have decreased and the formal 
quality of the fit, as measured with the x 2 and Q-parameter, 
is actually reduced. 

The scatter might be reduced if more accurate K-band 
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Figure 7. Deviations of our objects with respect to the baryonic 
mass vs. Vasymp TF relation, as a function of the colour excess of 
each galaxy with respect to the linear fit shown in figure[5] Symbol 
colours are the same as in figure [2] The dotted line shows the 
beha viour predicted by th e stellar population synthesis models 
from lBell fc de Jond l|200ll l. 



mass-to-light ratios become available for individual galaxies 
instead of the constant value of M*/Lk = 0.8 we assumed 
here. We investigate this possibility in figure where we 
plot the deviation of each galaxy with respect to the bary- 
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excess with respect to the linear fit shown in figure [5] As 
discussed in section l4.1.2l the global colour- magnitude trend 
only leads to a change in slope of the TF relation. The devi- 
ations from this global trend, however, are expected to give 
rise to additional scatter in the TF relation, hence the use 
of the colour excess on the abscissa. 

It is clear from figure [7] that there is little relation be- 
tween the offsets from the baryonic TF relation and the 
colours. In particular, the data points do not follow the trend 
as predicted on the basis of the ste llar populations synthesis 
models from lBell fc de Jond ([2001 ). Given the uncertainties 
in these models and in our own data, we cannot yet rule 
out the possibility that part of the scatter in the baryonic 
Tully-Fisher relation is related to variations in the mass-to- 
light ratios, but the current data suggest that this is not the 
main contribution. Thus, most of the scatter in the baryonic 
Tully-Fisher relation must be related to other observational 
effects, or be intrinsic. 

Note that the value of the \ 2 parameter for the baryonic 
Vksymp TF relation (1.92 or 1.72, depending on whether or 
not the outlying UGCs 3993 and 6787 are included in the fit) 
implies that roughly 50% of the scatter in this relation re- 
sults from the observational uncertainties in the data points. 
This leaves an upper limit of approximately 0.25 magnitude 
(« 25% luminosity) for the intrinsic TF scatter in our galax- 
ies. Given the fact that we have not included the uncertain- 
ties in the internal extinction corrections or the mass-to-light 
ratios in our errorbars and that there may well be additional 
uncertainties which we have ignored, it seems safe to assume 
that this is a conservative upper limit and that the intrinsic 
TF scatter is very small indeed. 



5 CONCLUDING REMARKS 

In this study, we have used K-band photometry and high- 
quality rotation curves to show that there is a 'break' 
in the Tully-Fisher relation around a rotation velocity of 
about 200 km s" 1 (equivalent to M C K w -23.75, or M* « 
0.7 x lO^M©), above which most galaxies rotate faster than 
expected (or equivalent, are less luminous). This kink is most 
pronounced in the traditional formulations of the TF re- 
lation, using stellar luminosities and the widths of the Hi 
profiles or the maximum rotational velocities, and we have 
shown that it is mainly due to massive early-type disk galax- 
ies with declining rotation curves, which lie systematically 
to the right of the relation defined by less massive and later- 
type spiral galaxies. The kink is also present, albeit to a 
lesser extent, in the relation with the asymptotic rotation 
velocity at large radii. In this case, there is no systematic 
offset between early- and late-type spirals, and the change 
in slope appears independent of morphological type. 

It is important to note that the 'elliptical galaxy Tully- 
Fisher relation' (using circular velocities derived from dy- 
namical models, rather than direct observations) shows a 
simil ar change in slope, at approximately the same loca- 
tion (|Gerhard et all 1200 it Ide Riicke et al.ll2007t ). Thus, the 
change in slope of the stellar luminosity vs. circular veloc- 
ity relation appears to hold universally and along the entire 
Hubble sequence. Note, however, that due to the limited 
spatial extent of the dynamical tracers in elliptical galaxies, 



it is not generally possible to measure the asymptotic rota- 
tion velocity in these systems. It can therefore not be ruled 
out that the kink in the elliptical galaxy TF relation would 
(partly) disappear if one could use the circular velocity at 
radii comparable to the outer edges of the Hi disks in spiral 
galaxies. 

A change in slope at the high-luminosity end has impor- 
tant consequences for the use of the Tully-Fisher relation as 
a tool for estimating distances to galaxies or for probing 
galaxy evolution. In recent years, many studies have ad- 
dressed the evolution of the Tully-Fisher relation on cos- 
mological timescales (i.e. between z « 1 and 0), although 
the results are still somewhat ambi guous. In their pioneer- 
ing studies, IVogt et"ail (1 19961 . Il997h reported that galaxies 
at redshift z ~ 1 were on average 0.6 mag brighter (B- 
band) than galaxies in the local un iverse. Other authors 
find much larger luminosity evolution l|Milvang- Jensen et al.l 
2003; iBamford et al.l |2006 | ). up to 1.5 - 2 magnitudes at 
redsh ifts of 0.25 - 0.45 |Rix et ail 1 19971 ; ISimard fc Pritchetl 
1998), whereas it has also been claimed recently that this 
is all due to observational effects, and that in reality, there 
has b een no luminosity evolution whatsoever l|Flores et al.l 
2006). Given the uncertainties in the measured zeropoint 
offsets, it is not surprising that me asurements of evolu - 
tion in the TF slo pe d isagree as well. Ziegler et al.l (2002), 
iBohm et all (|2004h and lBohm fc Zieglerl <|2006f > claimed that 
the evolution in the Tully-Fisher relation is luminosity de- 
pendent, with high mass galaxies (Vm ax > 150 km s -1 ) 
showing little or no evolution, but low mass galaxies be- 
ing up to 2 magnitu des brighter at high redshift. In contrast, 
IWeiner et all (120061 ) found that massive galaxies evolve more 
than fainter ones. 

Whatever concensus will be reached eventually, the 
presence of a change of slope in the local TF relation needs 
to be taken into account when interpreting the observed 
evolution. Our results indicate that high mass galaxies are 
under-luminous in the local universe, compared to a sim- 
ple extrapolation of the linear relation for lower-luminosity 
galaxies. This offset is strongest for early-type disks, i.e. 
galaxies with a high-surface brightness bulge, and for the 
TF relation based on the maximum rotation velocity Vmax- 
Crucially, high-redshift studies will often be biased towards 
such high-surface brightness galaxies. Moreover, when based 
on optical spectroscopic measurements, they will not be able 
to detect a declining rotation curve at large radii, and thus 



be forced to use Vm 



rather than V af 



These two se- 



lection effects, when applied to our own data, give rise to 
an offset of high-mass galaxies from the main relation, de- 
fined by low-mass galaxies with flat rotation curves, of 1 - 2 
magnitudes (middle panel in figure[2]). Thus, the luminosity 
evolution of high mass galaxies may be much larger than 
derived by the authors mentioned above. 

We have also shown that the 'kink' in the TF relation 
can be reduced in two different ways, namely 1) by using 
the asymptotic rotation velocity from the rotation curve as 
kinematic parameter and 2) by using the total baryonic mass 
(stars + gas) rather than luminosity. The first correction ap- 
pears to be the most important one for high mass galaxies 
and the TF relations using the asymptotic rotation veloc- 
ities show only a small deviation from linearity. However, 
the inclusion of the gas mass also improves the linearity and 
only when both refinements are used in conjunction does 
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the 'kink' disappear completely and is a linear relation re- 
covered. 

It is, again, interesting to note that the kink in the ellip- 
tical galaxy TF relation also disappears when using the total 
baryonic mass rather than stellar luminosity J Gerhard et al.l 
l200ll ; Ide Riicke et"aT]|2007l ) , so that galaxies along the entire 
Hubble sequence appear to follow the same, linear, baryonic 
TF relation. Our results seem a strong confirmation of the 
idea that the baryonic Tully-Fisher relation is fundamen- 
tally a relation between the mass of the dark matter haloes 
(which define V aS ymp) and the total baryonic mass in galaxies 
and that it holds universally, regardl ess of how the baryons 
are distributed within the haloes (cf. iMcGaueh et al 1 l200d : 
IVerheiienllioOll ). 
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